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Abstract 

In this work, we study the numerical optimization of nearest-neighbor concurrence of bipartite 
one and two dimensional lattices, as well as non bipartite two dimensional lattices. These systems 
are described in the framework of a tight-binding Hamiltonian while the optimization of concurrence 
was performed using genetic algorithms. Our results show that the concurrence of the optimized 
lattice structures is considerably higher than that of non optimized systems. In the case of one 
dimensional chains the concurrence is maximized when the system begins to dimerize, i.e. it 
undergoes a structural phase transition (Peierls distortion). This result is consistent with the 
idea that entanglement is maximal or shows a singularity near quantum phase transitions and 
that quantum entanglement cannot be freely shared between many objects (monogamy property). 
Moreover, the optimization of concurrence in two-dimensional bipartite and non bipartite lattices 
is achieved when the structures break into smaller subsystems, which are arranged in geometrically 
distinguishable configurations. This behavior is again related to the monogamy property. 

PACS numbers: 03.67.-a,03.65.Ud,73.43.Nq,71.10.Fd 



1 



I. INTRODUCTION 



Quantum entanglement is one of the most distinctive features in quantum mechanics 
and yet its properties are still not fully understood. This quantum resource is considered 
a key element for several quantum information and quantum computation proposals such 
as quantum teleportationi, superdense coding^, certain kinds of quantum key distribution 
schemes and quantum secret sharing protocols^i^. 

Recently, much research has been focused on a better understanding of quantum cor- 
relations in multiparticle systems^i^i^i^. A characteristic property that sets apart quantum 
correlations (or entanglement) from the classical ones is that entanglement cannot be freely 
shared among many objects. For the special case of three qubits A, B and C, this has 
been shown by Coffman et al.- using concurrence, a measure of entanglement taking values 
between and 1—. Coffman et al. showed that the sum of the squared concurrences AB 
and BC cannot be greater than unity, meaning that the degree of entanglement between 
A and B limits the entanglement between A and C . This property is called monogamy of 
entanglement. For an infinite chain of qubits sharing uniform entanglement with their first 
neighbors, Wootters^ found an upper limit for concurrence of 1/v^ (this limit, however, has 
not been yet proven for physical systems). Under certain conditions, Wootters reports a 
maximum concurrence of Cmax = 0.434467, a result which has also been confirmed in the 
case of rings with qubits when oo^. Diir et al.- have shown that in the case of 

three qubits, average concurrence between pairs of qubits of the \W) state are all equal to 
2/3, while Koashi et al.- have reported a maximum average concurrence between all pairs 
of qubits oi 2/N when all qubits in the system except one are in the state |0). 

On the other hand, it has been conjectured by several authors that entanglement can play 
an important role in quantum phase transitions (QPT)ii'i^ii^ii^ii^ii^ii^, which take place at 
absolute zero temperature^^. For example, Osborne et al.— have stated that in QPT, long- 
range correlations take place due to entanglement. In addition, Osterloh et al.— analyzed 
the behavior of entanglement near the critical point of the spin 1/2 model XY in a trans- 
verse magnetic field and found that in the region close to a quantum phase transition the 
derivatives of entanglement of formation (concurrence) obey a scaling law. Thus, there exists 
an intimate connection between entanglement, scaling and universality. Further, quantum 
entanglement has been used not only near quantum phase transitions to characterize them, 
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but also to obtain a better description of experimental measurements of specific heat and 
magnetic susceptibility of dilute solution of Ising dipoles^. Additionally, we have recently 
shown that, in the case of rings with off-diagonal disorder, an increase in disorder strength 
results in enhancement of the nearest-neighbor (NN) concurrence with respect to the per- 
fectly ordered ring2S. This enhancement of concurrence hints that quantum entanglement 
might be a good indicator of anomalies in the wave functions and density of states in systems 
where no quantum phase transitions take place. 

In the present paper, we report maximization of the concurrence using a computational 
optimization method. The calculations were done on the basis of an electronic model Hamil- 
tonian which contains parameters depending on the lattice structure. Such parameters were 
optimized in order to maximize the concurrence. This means that, in general, the systems 
considered will not be translationally invariant after optimization, since the structure can 
be dramatically changed by the genetic algorithm. We present, specifically, calculations 
on systems described by a tight-binding model in one and two dimensional lattices. Our 
results can serve as a basis for the design of structures which maximize nearest-neighbor 
entanglement. 

This paper is organized in the following manner: in section II the Hamiltonian is pre- 
sented, as well as the most important formulas used to quantify concurrence and the key 
steps in the optimization procedure. Results for one and two dimensional bipartite lattices, 
as well as non-bipartite two dimensional lattices will be presented in section III. Finally, a 
summary of our results can be found in section IV. 

II. THEORY 

We consider electronic systems described by a tight-binding Hamiltonian of the form 



For simplicity, we consider spinless electrons. In Eq. ([T]) cf (cj) is the usual creation (anni- 
hilation) operator of a spinless electron at site i and tij is the hopping integral between NN 
sites i and j. 

It should be noted that the magnitude of the hopping elements t^j can be easily related 
to the lattice structure. A large tij indicates a small interatomic distance, which causes a 
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large overlap between the wave functions localized on the sites i and j. Analogously, a small 
tij represents a large interatomic distance, tij = means that sites i and j are not nearest 
neighbors. 

To calculate the entanglement of formation between pairs of sites we employed the concept 
of concurrence^ and proceeded in a similar way as Zanardi and coworkers did for the case 
of translationally invariant chains^i. 

For a system described by the Hamiltonian of Eq. ([1]), the concurrence between sites i 
and j has the forn>^ 

Q,- = 2max{0,|^| -y^}, (2) 

where 



V = I - (hi) - {hj) + {hihj) (3) 
y = (nirij) (4) 
z={c^c,). (5) 

Note that i and j do not need to be nearest neighbors. 

To maximize NN concurrence, we use genetic algorithms (GAs), a technique developed 
by J. H. Holland and his students at the University of Michigan in the 1960s and 1970s. 
Holland's goal was not to design algorithms to solve specific problems, but rather to formally 
study the phenomenon of adaptation as it occurs in nature and to develop ways in which 
the mechanisms of natural adaptation can be used for numerical optimization. 

In this optimization technique, the characters or numbers representing the solution to a 
problem are stored in a string called "chromosome" . Nature-inspired operators of crossover 
and mutation are applied to a population of chromosomes and those individuals that repre- 
sent the best solution to the problem are given more probability of being chosen for the next 
generation. By iterating this process, it is possible to obtain a very good solution without 
having to explore the entire solution space. 

For the particular problem treated in this paper, we propose first families of hopping 
integrals {tij}, which are then fine-tuned by means of the genetic algorithm. Our goal is 
to obtain the hopping elements from the Hamiltonian such that its ground state function 
represents the maximum average NN concurrence. 
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We have used genetic algorithms in a previous paper—, where we have optimized the total 
concurrence (Ctotai) of one dimensional chains and two dimensional lattices described by a 
tight-binding Hamiltonian with open and periodic boundary conditions. Ctotai was taken as 
an average over all concurrences between every site i with the — 1 remaining sites. 

In this paper, however, we focus only on optimization of first-neighbor concurrence, which 
can be calculated by means of the following expression: 



where z is the number of nearest neighbors. Cnn is our fitness function for the application 
of the genetic algorithm. 

It is important to point out that a system with translational invariance is described by all 
the tij elements from Eq.[T] having the same value, for example, tij = — 1. This configuration 
leads to concurrence values that are the same for each pair of nearest neighbors, a case which 
has been previously studied by several authors^'^'^. 

The genetic pseudo-algorithm employed in this work has been given in detail elsewhere^^. 
Briefly, it consists in the following steps: 

1. We consider each tij of the Hamiltonian matrix as a gene and the array of these genes 

ClS db chromosome. 

2. Allocate two arrays, "generationO" and "generation!" composed of chromosomes. 

3. Allocate a chromosome "best" with fitness 0.0. 

4. For a given band filling, we repeat the following steps: 



• Initialize "generationO" with random values in the range (—5,0). 

• For a given number of generations repeat: 

— Decode each chromosome in "generationO" into a Hamiltonian matrix, diag- 
onalize it and calculate the average concurrence between only nearest neigh- 
bors of the system using Eq. [6l In other words, we calculate the fitness Cnn 
for each individual in "generationO" . 

— Compare "best" with the fittest individual in "generationO" and substitute 
"best" if the latter has larger fitness. 




(6) 



i=i j=i 
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— Choose chromosomes with a probabihty proportional to its fitness (selection 
operator) and copy it to "generationl". 

— Use crossover and mutation operators on chromosomes in "generationl" to 
create new chromosomes. 

— Make "generationO" equal to "generationl". 

• Print "best" in an output file. 

III. RESULTS AND DISCUSSION 

In this section, we present our results of NN concurrence optimization using GAs. The 
effect of bipartite and non-bipartite systems with periodic boundary conditions is addressed. 
We recall that a lattice is considered bipartite if it can be separated in two sublattices A 
and B such that all first neighbors of the sites of sublattice A are sites from sublattice B 
and vice versa. Biparticity or non-biparticity has important consequences in the physical 
properties of a lattice. For example, the concurrence of a bipartite lattice is symmetric 
around half band filling while non-biparticity is responsible for magnetic frustration in spin 
systems, that is, the impossibility of minimizing energy for each pair of spins in the lattice. 

A. Bipartite Systems 

First, we focus on finite one dimensional chains of up to = 50 sites with periodic 
boundary conditions. Results are shown in Fig. [H We considered 500 chromosomes (i.e. 
a collection of hopping integrals from the hamiltonian matrix) in each population, with 
starting values randomly chosen in the range (—5, 0). The algorithm was left to evolve 4000 
generations in order to obtain the individuals with the best fitness value. This procedure 
was undertaken for each band filling (x = n/N). From the figure it can be seen that there 
exists symmetry around half band filling, indicating that it is possible to focus only on the 
< a: < 0.5 region when analyzing concurrence in bipartite lattices. 

Optimized concurrence shows two distinctive features: (a) the growth of concurrence as 
function of band filling is quasi-linear for x < 0.25. This result is due to the fact that in lower 
band fillings, the probability of electronic collision is low, and the system is able to optimize 
concurrence without considering Pauli repulsion, (b) For x > 0.25 the growth of Cnn as 
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function of x is quasi-parabolic, and electronic exchange due to the Pauli principle comes 
into play. Comparison between results obtained using GAs and the ordered case -with all 
tij taken equal- shows that for low band filling (x < 0.1), both concurrence curves display 
a similar behavior and the impact of optimization starts to be noticeable for x > 0.1, with 
a maximum increase in concurrence with respect to the ordered case of 47% (0.5 vs 0.3392) 
at half band filling. 

In order to understand the behavior of optimized concurrence, we have analyzed those 
chromosomes with best fitness values for each x. Results show that, as band filling x is 
increased, the ring starts to break into smaller chains (i.e., some hopping elements become 
tij ^ 0). For example, while for n = 1 all tij are similar and close to —4, in the next filling, 
n = 2, two noncontiguous elements tij differ from the rest and are now close to 0, suggesting a 
separation into two subchains. This division increases with x until half band filling (x = 0.5), 
a point at which the system alternates a short bond {tij ~ — 5) with a long one (tij ~ 0) 
as can be seen depicted in Fig. O In other words, to increase the average concurrence, the 
system begins to dimerize, i.e, the chain undergoes a structural phase transition due to a 
Peierls instability^'^^'^. This is a metal-insulator transition occurring in one-dimensional 
metals, where the doubling of the unit cell leads to a decrease in the kinetic energy of the 
system. Note that dimerization of the chain was obtained in a natural manner, since all 
initial matrix elements were initialized by random values. In fact, the maximum concurrence 
at half band filling corresponds to a system which has completely dimerized. This shows a 
possible connection between the increase in concurrence and the Peierls instability. In order 
to study this connection, we analyzed the concurrence and their derivatives for a dimerized 
chain. In this case, we consider that the nearest-neighbor hopping integrals tij of Eq. [T]take 
the values t2n,2n.+i = 1+a and t2n-i,2n. = 1— a, where a G [0, 1] is the dimerization parameter. 
This one-particle dimerized Hamiltonian can be diagonalized by using the transformations 
proposed by Su, Schrieffer and Heeger— . The goal of analyzing the concurrence and its 
derivative is based on the fact that a discontinuity (singularity) in the (derivative) ground- 
state concurrence has been associated to a first (second) order QPT— >^>^>^>^. However, 
it has been shown that the relation between QPT and non-analyticity in the concurrence 
is not one-to-oneiS. A one-to-one connection can be assumed though, when we consider 
QPTs characterized by nonanalytic behavior in the derivatives of the ground-sate energy 
and we exclude artificial and accidental occurrences of non-analyticities in the ground-state 
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concurrence and its derivative^. In Fig. [3l we present the ground-state concurrence of a 
dimerized chain, C2n,2n+i and C2n-i,2n,, and its derivative, C2„2n+i ^'2n-i2n^ where C 
means dC/da, as a function of a. Note that C2n,2n+i = 2.0*max{0, 72n,2n+i + 72n 2n+i "0-25} 
and C2n-i,2n = 2.0 * max{0,72n-i,2n + 7L-i,2n - 0.25>2a, where 72n,2n+i = (4nC2n+i) and 
72n-i,2n = (cln-i'^'in) are the one-particle density-matrix elements or bond orders between 
NN and can be calculated analytically^. These bond orders are continuous functions of a, 
the first one ranging from 72n,2n+i = 0.318310 (a = 0) to 0.5 {a = 1.0) and the second one 
from 72n-i,2n = 0.318310 (a = 0) to 0.0. (a = 1.0). Therefore, the discontinuity obtained 
for C2n-i,2n at « ~ 0.138 is not related to a critical point. Clearly, this discontinuity is 
artificial and comes from the particular definition of the concurrence in Eq. 

In Fig. 3 one observes that 2n+i C2„_i 2n present singularities at the limit a — *• 
ac = 0, where the Peierls instability occurs, which should be related with a second order 
QPT—. In order to investigate if the singularity is related with a second order QPT, we 
write the ground state energy (and their derivatives) as a function of the NN density- 
matrix elements (and their derivatives). The ground-state energy per site is given by -Egg = 
-(1 + a)72„,2n+i - (1 - a)72n-i,2n, its first derivative by dE^^/da = -{1 + a)d^2n,2n+i/ da - 
(1 - a)d'^2n-i,2n/da - 72n,2n+i + 72n-i,2n and its sccoud derivative by d'^E^^/d'^a = -(1 + 
a)d^-i2n,2n+i/d^a - (1 - a)d^-i2n-i,2n/d^a - 2rf72n,2n+i/c?a + 2d-i2n-i,2n/da. These second 
derivatives presents a singularity at Uc-, which is a manifestation of a second order QPT. 

On the other hand, the maximal value of the Cnn is obtained when the chain is completely 
dimerized (a = 1), i.e. the chain has been transformed into N/2 singlet Bell states. This 
complete dimerization can be related to the ideas of monogamy between three systems^: the 
optimized system prefers to have maximum concurrence (C=l) between subsystems A and 
B while lowering that between B and C (C=0). 

The effect of dimensionality over Cnn is shown in Fig. HI where calculations for a square 
lattice 6x6 using a population of 700 individuals were performed. We considered 10,000 
generations for each band filling. From the figure, it is possible to observe again two clearly 
remarkable traits for two ranges of x, as in the one dimensional low band filling 

range (x < 0.25) where the increase of Cnn as function of x is quasi-linear due to low 
collision probability between particles, and a second range x > 0.25 where the Pauli exclusion 
principle starts to play an important role. Cnn as function of x has a quasi-parabolic behavior 
in the latter range, decreasing first in the range 0.25 < x < 0.375 and then increasing 
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for 0.375 < X < 0.5. It is also noteworthy that both regions have the same maxima of 
Cnn ~ 0.25, one at x = 0.25 and the other at x = 0.5. Comparison of Cnn between the 
optimized case and an ordered, large lattice shows that the greatest difference occurs at 
half band filling. Moreover, we observe that maxima Cnn obtained for the square lattice is 
smaller than the maximum Cnn for the ID ring. 

A study of the fittest chromosomes for those band fillings where Cnn is a maximum shows 
the following behavior: for x = 0.25, the short bonds {tij ~ —5) form nine squares which are 



separated by long bonds {tij ^ 0) (Fig. 5(a) ). In each square, a spinless electron can be found. 
This result is related to the definition of generalized \ Wn) states, \Wn) = {1/ VN)\N — 1, 1) 
where the state — 1, 1) denote the totally symmetric state with — 1 zeros and one 1. 
It has been shown that these \Wn) states are very robust against particle losses and that 
the concurrence of two qubits can be determined to be Cij = 2/N—. In the case of = 4, 
we have that 11^4) = (1/^4) [|0001) + |0010) + |0100) + |1000)] and Cij=0.5 for two NN 
sites of a square. Using this Cij value for two NN sites of a square, it is easy to show that 
Cnn=0.25 (see Eq. [6]). On the other hand, for x = 0.5 the system tries once again to form 



dimerized states, separated with long bonds (Fig. 5(b)). The decrease of Cnn with respect 
to the one dimensional ring is related to the difference in the number of neighbors between 
both systems. In the one dimensional systems, Cnn is obtained by averaging between two 
nearest-neighbors whereas in the square lattice this average requires four nearest-neighbors. 



B. Non Bipartite Systems 

In order to study the effect of non-biparticity on the concurrence, we have considered three 
non-bipartite lattices: the Kagome lattice, the maple leaf or Betts lattice and the triangular 
lattice. Note that the triangular lattice is less bipartite than the two other ones due to its 
larger number of triangular bonds. This has important consequences in antiferromagnetic 
systems^^. 

The basic common feature of the results presented so far is that successful individuals in 
the linear chain and square lattices form open and tightly closed binary subsystems. Thus, 
we exploited this fact in order to speed up the procedure and increase the accuracy of the GA 
calculations. We initialized a hundred individuals of the population to 0.0 and a hundred to 
—5.0. The motivation for performing this step is that this short and long bonds will spread 
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to other individuals through the crossover operator, and that this kind of bonds will help 
increase its concurrence. 

We now apply our calculations to a periodic Kagome lattice of 48 sites with periodic 
boundary conditions. Results are shown in Fig. [Hi We considered 1200 individuals in the 
population, which was left to evolve over 3000 generations. Quasi-linear behavior of Cnn 
with X can be seen in regions with low electronic or hole density, namely x < 0.2 and x > 0.8. 
This behavior corresponds to low probability of collision, as has already been discussed, and 
is independent from the kind of lattice or whether it presents triangular bonds or not. 
The first effect of non-biparticity is reflected in the lack of symmetry of Cnn around half 
band filling. Note that Fig. E] presents only one distinguishable maximum located -contrary 
to the case bipartite lattices-, out of half band filling, at x ^ 1/3. Comparison with 
respect to the ordered case shows the relevance of the optimization technique, particularly 
for 0.6 < X < 0.8, where the initially ordered structure displays a very low concurrence. 
Moreover, there is a coincidence in the localization of the maximum for both the optimized 
and ordered lattices. The best individual corresponding to the maximum of the numerical 



optimization approach can be seen in Fig. 7(a) , where one can observe that a geometrical 
configuration favoring triangular loops has been formed. It is clear again that this trend to 
form triangular loops is related to W states, \W) = (1/^3) [|100) + |010) + |001)], where 
Cij = 2/3 for two NN sites forming the triangle. To further extend the optimized results. 



we now introduce a design as seen in Fig. 7(b) A new concurrence curve for the Kagome 
lattice is presented in Fig. El In this case, besides the hundred individuals initialized to 0.0 



and —5.0 each, all other individuals were initialized with the proposed design of Fig. 7(b) 



As can be seen, the concurrence is not only increased at the maximum, but the proposal 
improved the overall behavior of the curve as well. 

Considering the hint given by the analysis of the Kagome lattice, we ran the optimization 
algorithm for the maple-leaf or Betts lattice in order to find an optimum structure design. 
Fig. [8] shows the best individual of an initial run of the algorithm corresponding to the 
maximum at x = 0.33 [Fig. 8(a)| , as well as a proposal based on this individual [Fig. 8(b)| . 
Fig. [9] shows the results for the ordered CclSG, clS well as both runs of the algorithm. As 
expected, the ordered case performed worse than any of the optimized cases. Once again. 



the run based on proposal shown in Fig. 8(b) performed better than the random-based case. 
Finally, in the case of the triangular lattice, the typical behavior found in the other non- 
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bipartite lattices can be readily observed: a single maximum located near x = 1/3 and a 
much better performance with respect with the ordered case (see Fig. [TOj) . As in the other 
cases, we then proceeded to develop new proposed configurations by looking at successful 
individuals and being inspired by W states. Two new designs were chosen, seen in Figures 



11(a) and 11(b) 



The best curve for Cnn can be seen in Fig. [TOl and is basically a combination of the dif- 
ferent curves obtained when initializing with each proposal. Although there is a fair increase 
in the overall maximum (Cnn = 0.2222), it is still lower than that of the Kagome lattice 
(Cnn = 0.3333) and maple leaf lattice (Cnn = 0.2666). This is because in the triangular 
lattice there are more bonds, and consequently each site must share its entanglement with 
more neighbors. 



IV. SUMMARY AND CONCLUSIONS 



In this work, we have used OAs to maximize the nearest-neighbor average concurrence 
of systems by tuning the nearest-neighbor hopping integrals of a tight binding Hamiltonian. 
The optimization of entanglement has been performed for one- and two dimensional bipartite 
systems as well as for two dimensional non bipartite systems. The results show that the 
concurrence of the optimized systems is very large in comparison with the ordered structures. 
This increase in the concurrence is understood and interpreted by analyzing the optimized 
nearest-neighbor hopping integrals. In general, we found certain tendencies of periodical 
systems to break into smaller subsystems. This is achieved in a natural manner by the 
system by making the hopping integrals evolve in such a way that the absolute value of 
some integrals is high in some cases {\tij\ ~ 5) and very small in others {tij ~ 0). 

This results are related to the fact that quantum entanglement, in contrast with classical 
correlations, cannot be freely shared between many objects. This quantum correlations 
property — monogamy — , is clearly noticeable in the case of the periodic ring at x = 0.5, 
where the n-th site is completely entangled with the site n + 1 (Cn,n+i — 1) while on the 
other hand it is completely unentangled with the site n — 1 (Cn,n-i — 0). Moreover, results 
at X = 0.5 shows that, in order to maximize concurrence, the system undergoes a structural 
transition, the Peierls distortion. 

Finally, it is worth mentioning that the relationship between the ID tight binding Hamil- 
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tonian [Eq. ([T])] and the XX chain of 1/2 spin given by the following expression 

Na J 

j 

where Jj is the coupling constant. This relationship can be shown by using the Jordan- 
Wigner transformation^^ and therefore, it can be inferred that our results for the ID tight 
binding systems can be used to describe the maximization of the Cnn concurrence for systems 
modeled by the XX chain of 1/2 spin. 
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One dimensional ring 

N = 50; gen. = 4000; pop. = 500 




FIG. 1: (Color online) Nearest-neighbor concurrence Cnn of chains with periodic boundary con- 
ditions as a function of band filling x = n/N. 
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FIG. 2: (Color online) Graphical representation of a one-dimensional ring with 32 sites in the 
lattice structure yielding the maximal nearest-neighbor concurrence (see text). Parameters: 5000 
generations; population=600. Half band filling. Three cases are distinguished by color: light grey 
{\tij\ < 1), grey (1 < < 3) and black (3 < < 5). 
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FIG. 3: (Color online) (a) Nearest-neighbor ground-state concurrences Cj^j+i of the dimerized chain 
(b) and their derivatives as a function of the dimerization parameter a. 
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Square lattice 

N =36; gen. = 10,000; pop. = 700 




X = n/N 

FIG. 4: (Color online) Nearest-neighbor concurrence Cnn of a 6 x 6 square lattice as a function of 
band filling x. The solid line refers to the optimized system. The dashed line shows the concurrence 
of a large perfectly periodic square lattice {N = 1600). 




(a)Band filling = (b)Band filling = 

25% 50% 

FIG. 5: Graphical representation of the optimized structures obtained starting from a square 
lattice of 36 (6 x 6) sites at (a) a quarter band filling and (b) half band filling. Parameters: 10000 
generations; population=700. Three cases are distinguished by color: light grey {\tij\ < 1), grey 
(1 < \tij\ < 3) and black (3 < \tij\ < 5). 
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Kagome lattice 

N = 48; gen. = 3,000; pop. = 1 ,200 




0.4 0.6 0.£ 

X = n/N 



FIG. 6: (Color online) Nearest-neighbor concurrence Cnn of small periodic Kagome lattice as a 
function of band filling x. In addition, nearest neighbor concurrence Cnn for the Kagome lattice 
were calculated using the proposed optimal structure of Fig. |7(b)| 
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(b)Proposed optimal 
structure 



FIG. 7: Graphical representation of a Kagome lattice of 48 sites, (a) represents the best individual 
using genetic algorithms at x = 33%. (b) depicts a proposed optimal structure based on (a). 
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(a)Best individual for Betts lattice (b)Proposed optimal structure 

FIG. 8: Graphical representation of a Betts lattice of 48 sites, (a) represents the best individual 
using genetic algorithms at x = 33%. (b) depicts a proposed optimal structure based on (a). 
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Belts lattice 

N = 54; gen. = 3,000; pop. = 1 ,200 
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FIG. 9: (Color online) Nearest-neighbor concurrence Cnn for the Betts lattice as a function of 
X. The figure shows Cnn for the ordered case, the optimized case using genetic algorithms and, 
finally, the case where the generations were initialized at each band filling using the proposed 
optimal structure of Fig. |8(b)[ 



Triangular lattice 

N = 36; gen. = 3,000; pop. = 1 ,500 
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FIG. 10: (Color online) Nearest-neighbor concurrence Cnn for the triangular lattice as a function 
of X. The nearest-neighbor concurrence Cnn was calculated for the ordered case, the optimized 
case using genetic algorithms, the proposed optimal structures and the case combining results from 
proposals. 
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I3^6^jc^ ^^^^^ 

(a)Proposcd optimal (b)Proposcd optimal 

structure "Triangle" structure "Diamond" 

FIG. 11: Proposed optimal structures based on the best individual obtained using genetic algo- 
rithms. The "triangle" structure is based on the best individual at x = 33% whereas the diamond 
is inpired on the best individual at x = 0.69. 
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